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A new regime of fast magnetic reconnection with an out-of-plane (guide) magnetic field is reported 
in which the key role is played by an electron pressure anisotropy described by the Chew-Goldberger- 
Low gyrotropic equations of state in the generalized Ohm’s law, which even dominates the Hall term. 
A description of the physical cause of this behavior is provided and two-dimensional fluid simulations 
are used to confirm the results. The electron pressure anisotropy causes the out-of-plane magnetic 
field to develop a quadrupole structure of opposite polarity to the Hall magnetic field and gives rise 
to dispersive waves. In addition to being important for understanding what causes reconnection 
to be fast, this mechanism should dominate in plasmas with low plasma beta and a high in-plane 
plasma beta with electron temperature comparable to or larger than ion temperature, so it could 
be relevant in the solar wind and some tokamaks. 


Magnetic reconnection allows for large-scale conver¬ 
sion of magnetic energy into kinetic energy and heat by 
changing magnetic topology. It occurs in a wide range of 
environments, such as solar eruptions, planetary magne¬ 
tospheres, fusion devices, and astrophysical settings [ 1 ]. 
One key unsolved problem is what determines the rate 
that reconnection proceeds mig. 

In simplified two-dimensional (2D) systems often em¬ 
ployed in simulations, the reconnection rate E is deter¬ 
mined by the aspect ratio of the current sheet, but it is 
not understood what controls its length. In collisional 
plasmas, current layers are elongated mu which make 
collisional reconnection relatively slow. For less colli¬ 
sional 2D systems, elongated layers break and produce 
secondary islands ism, giving normalized reconnection 
rates of O.OI jOUTT]. However, this is ten times slower 
than the fastest rates seen in simulations nMz]. 

The GEM Challenge [18] showed that the Hall term, 
when active, is sufficient to produce short current layers 
with E 0.1. The interpretation of this is still under 
debate [TM25] . One can ask - do other mechanisms limit 
the length of current layers that could help explain what 
causes fast reconnection? 

In this Letter, we report that fast reconnection can be 
caused by electron pressure anisotropy using the Chew- 
Goldberger-Low (CGL) equations of state [26] in the gen¬ 
eralized Ohm’s law. This has not been seen previously 
because (a) most simulations use no out-of-plane (guide) 
magnetic field, but this mechanism requires one and (b) 
previous fluid simulations included pressure anisotropies 
only in the momentum equation, which does not pro¬ 
duce fast reconnection [27H29]. This result is distin¬ 
guished from known results that off-diagonal elements 
of the electron pressure tensor balance the reconnection 
electric field at the reconnection site [3QH32] and agy- 
rotropies contribute near the reconnection site [SS] [34] . 
Neither effect is present for the GGL equations because 
the pressure tensor is gyrotropic. As with the Hall effect, 
gyrotropic pressure does not break the frozen-in condi¬ 


tion [35]. Nonetheless, it plays a crucial role in allowing 
fast reconnection in this regime. 

Gyrotropic pressures are different parallel p|| and per¬ 
pendicular p± to the magnetic field [36]. The GGL, or 
double adiabatic, equations of state [26] follow rigorously 
from kinetic theory in the ideal limit (no heat conduc¬ 
tion) with strong enough magnetic fields so particles are 
magnetized. Previous studies treated gyrotropic pres¬ 
sures in tearing instabilities [37H42] . Electron pressure 
anisotropies have garnered interest lately since they are 
self-generated by reconnecting magnetic field lines [43] . 
The resulting equations of state [44] are valid for guide 
fields no stronger than the reconnecting magnetic field. 

Simulations are carried out using the two-fluid code 
E3D [45] modified to include gyrotropic pressures. The 
code updates the continuity equation, momentum equa¬ 
tion, Earaday’s law, and pressure equations. The electric 
field E is given by the generalized Ohm’s law. 
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where v is velocity, B is magnetic field, J is current den¬ 
sity, n is number density, e is proton charge, Pe is the 
electron pressure tensor, 77 is resistivity, rrie is electron 
mass, and each term on the right can be turned off, in¬ 
cluding the Hall term J x B/nec. The momentum equa¬ 
tion is 
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( 2 ) 


where e = I — 47 r(p|| —p±)lB‘^, I is the identity tensor, p 
is mass density, and p is total (electron plus ion) pressure. 

When pressure anisotropies are used, we employ the 
GGL equations, equivalent to j and Pa^/pB be¬ 

ing constants [26] for species a. We write them as evolu¬ 
tion equations {see Eqs. (17) and (18) in Ref. [l 6 ]} with 
E • J omitted for simplicity. The numerical implementa¬ 
tion is benchmarked using Alfven waves and the firehose 
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and mirror instabilities. For isotropic plasmas, Palp^^^ 
= constant. 

All quantities are normalized: magnetic fields to the 
reconnecting magnetic field densities to the value po 
far from the current sheet, velocities to the Alfven speed 
ca = BqI (47rpo)^^^5 lengths to the ion inertial length di = 
cjujpi, electric fields to caBqIc^ resistivities to dTrc^d^/c^, 
and pressures to Bq/A n. The simulation size is Lx xLy = 
204.8 X 102.4 in a doubly periodic domain with 4096 x 
2048 cells. This system is large enough that boundaries 
do not play a role; a steady state prevails for an extended 
time. 

The initial reconnecting magnetic field is Bx{y) = 
tanh[(^ + I/^/4)/0.5] — tanh[(^ — I/^/4)/0.5] — 1. Un¬ 
less otherwise stated, the guide field is large at 5 and 
increases at the current sheet to balance total pressure. 
The density p = 1 and pressure Pcr = 5 are initially uni¬ 
form (pcr± = -Peril = 5 for anisotropic). When electron 
pressure is evolved, ions are cold, and vice versa. 

All simulations use rrie = mi/2b unless otherwise 
stated, which is acceptable because E is insensitive to 
me mi and length scales for the ions {cs/Lld — 0.7) 
and electrons {de = c/ujpe = 0.2) are sufficiently sep¬ 
arated [48]. The resistivity is 0.005, chosen so that if 
reconnection is Sweet-Parker-like, the layer thickness is 
(pI/a^/4)^/^ 0.5 which makes it marginal against sec¬ 

ondary islands [7]. Reconnection is seeded using a co¬ 
herent magnetic perturbation of amplitude 0.014. Initial 
random velocity perturbations of amplitude 0.04 break 
symmetry. The equations employ fourth order diffusion 
with coefficient = 2.5 x 10“^ to damp noise at the 
grid. 

Benchmark simulations using two-fluid simulations 
(with the Hall term, electron inertia, and isotropic elec¬ 
tron pressure) reveal a well-known open exhaust, as 
shown by the out-of-plane current density Jz in Fig.j^a). 
Various simulations are then performed without the Hall 
term. When the CGL equations are used on the electrons 
(which we call eCGL), an open exhaust occurs (panel 
b). Panel (c) is for the same system but with rrie = 0, 
showing that electron inertia does not cause the open ex¬ 
haust. Panel (d) is when the CGL equations are used on 
the ions (which we call iCGL). The current sheet is elon¬ 
gated like in Sweet-Parker reconnection. To further iden¬ 
tify the key physics, a simulation of an unphysical system 
is tested: the electron pressure anisotropy is included in 
the momentum equation [Eq. Q] but not in generalized 
Ohm’s law [Eq. Q]. The result is an elongated current 
sheet (panel e). The three with open exhausts are fast, 
E 0.06 — 0.1, while the elongated sheets give the Sweet- 
Parker rate of 0.01. The thickness of the current sheets 
in (a) and (b) are near 0.2, showing that layers in eCGL 
go down to de as in two-fluid reconnection. In contrast, 
the layer thickness for (d) and (e) is 0.525 and 0.6 (the 
Sweet-Parker thickness). The conclusion is clear: the 
eCGL equations give rise to fast reconnection even with 
no Hall term, and the key physics is the electron pressure 
gradient in generalized Ohm’s law. 

The physics when electron pressure anisotropies dom¬ 
inate bears similarities to Hall reconnection with a guide 
field |49|. The 2: component of Eq. 0 in terms of the 
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FIG. 1: (Color) Out-of-plane current density Jz using various 
models: (a) two-fluid, (b) eCGL, (c) eCGL with no electron 
inertia, (d) iCGL, (e) eCGL in the momentum equation only. 


flux function defined as Ez = (\/c)d^/dt is 


- T) •VV’ = -4:(V-Pe).+ + V>. 
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ne 


ne 
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(3) 

With the Hall term present, the left side reveals mag¬ 
netic flux is convected by electrons, so electrons carrying 
the current drag the reconnecting magnetic field out of 
the plane [50]. This produces a quadrupole out-of-plane 
magnetic field Bz m, shown for the two-fluid simulation 
in Fig. I^a). With a strong guide field, the gas pressure 
(not shown) develops a quadrupole with opposite polar¬ 
ity to maintain total pressure balance [49] . 

Without the Hall term, Eq. 0 implies magnetic flux 
is convected by ions [52]. The magnetic field is dragged 
out of the plane by ions, giving a Bz quadrupole with 
opposite polarity as in Hall reconnection, displayed for 
the eCGL simulation in Eig. [^b). (An instability is vis¬ 
ible in the exhaust. The system is not firehose or mirror 
unstable; it is likely a drift instability.) To balance total 
pressure, Pe± develops a quadrupole of opposite polarity, 
displayed in Eig. [^c). The density (not shown) develops 
a quadrupole like that of Pe±- This requires a parallel 
electric field pointing from low density to high, which 
comes from a parallel electron pressure with high Pe\\ in 
regions of low Pe± so Pe| | has a quadrupole of opposite po¬ 
larity as Pe-L, shown in Eig. [^d). Thus, an electron pres¬ 
sure anisotropy is self-generated. It contributes to the 
reconnection electric field as Epe^z = — (l/'ne)(V • Pe)z = 
— (l/ne)(B • V)[(pe± — Pe\\)Bz/B‘^]j plotted as the solid 
line in a vertical cut through the X-line in Eig. [^e). Eor 
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FIG. 2: (Color) (a) Out-of-plane magnetic field Bz in two- 
fluid reconnection. For eCGL simulation, (b) out-of-plane 
magnetic field Bz^ (c) perpendicular electron pressure Pe±, 
and (d) parallel electron pressure Pe||- Reconnection electric 
field contributions in the (e) eCGL, two-fluid, and iCGL sim¬ 
ulations and (f) a simulation with eCGL and the Hall term. 

comparison, the dashed line shows the Hall electric field 
Ehz in the two-fluid simulation and the dash-dot line 
shows the resistive electric field in the iCGL simu¬ 
lation. The structure of Epe^z is similar to the known 
Ehz profile in two-fluid reconnection. Note that eCGL 
with rUe = 0 also has quadrupoles, but iCGL with slow 


reconnection does not. 

The guide field is key to the physics. If it is too 
large, the ion Larmor radius falls below electron or resis¬ 
tive scales which prevents fast reconnection, analogous 
to Hall reconnection [53]. If it is too small, the pres¬ 
sure change due to the Bz quadrupole is small, so the 
effect in the previous paragraph is negligible. We quan¬ 
tify this by finding when Ep^^z is dominated by other 
contributions to Ohm’s law, which for the present sim¬ 
ulations is the resistive term. A scaling analysis gives 
Epe^zlEp ^ 0.1/3eo(5^/5rec)/(27^c^/47rcAG^i), where 0.1 
is E for fast reconnection, [3eo = SirpeoKB^^^ ^g) 
the electron plasma beta, and Brec and Bg are the re¬ 
connecting and guide field strengths. This ratio is small 
if Bg is sufficiently big or small. We confirm this in 
simulations with pe± = Pe\\ = I? Ihe predicted range is 
0.05 Bg 20. (Formally, the CGL model is invalid for 
small Bg^ so this tests fundamental physics independent 
of the appropriateness of the CGL model.) We find re¬ 
connection is Sweet-Parker-like for Bg = 0 and 0.1, has a 
short current layer with ^ 0.03 for a transitional guide 
field Bg = 0.25, is fast with E 0.05 for Bg = 0.5, 5, and 
7.5, and is again Sweet-Parker-like for Bg = 15. These 
results agree with the prediction. 

This system yields an interesting way to study the 
cause of fast reconnection. It was proposed that recon¬ 
nection is fast if linear perturbations to a homogeneous 
equilibrium permit dispersive waves with faster phase 
speeds at smaller scales m, such as the whistler or ki¬ 
netic Alfven wave in Hall-MHD. This has been contro¬ 
versial because reconnecting fields are not homogeneous. 

We present the linear theory of a plasma with pressure 
anisotropies, the Hall term, and electron inertia. Rather 
than using CGL, we generalize by taking ion and elec¬ 
tron pressures to be arbitrary functions of p and 5, z.e., 
= PcjAp^B) and p^\\ = p^\\{p,B). This captures 
adiabatic, CGL, and Egedal’s equations of state [43l 1^ . 
The dispersion relation relating the frequency uj to the 
wavevector k is 

- Co = 0, (4) 

where 
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Here, D = 1 -\- k‘^d‘1^ pb = H^/Stt is the magnetic 
pressure, e = 1 — 47r(p||o — p±o)/Bq is the equilib¬ 
rium anisotropy parameter for the total pressure, eg is 
similarly defined for the electrons, d{Ai, A 2 )/d{x,y) = 
{dAi/dx){dA 2 /dy) — {dAi/dy){dA 2 /dx) is a Poisson 
bracket-type operator, and the 0 subscript denotes equi¬ 
librium quantities. This reduces to known results in 
the limits of anisotropic-MHD with the CGL equations 
{da- O^PaxIpB = constant, p^\\B^ jp^ = constant) [55] 
and isotropic two-fluid (pcr|| =Pcr±) [54] . 

We find pressure anisotropies introduce dispersive 
waves even when the Hall term is absent. All terms 
in Eq. Q with df give dispersive waves. For the 
high cj, high k with e = eg = 1 limit, bP^ 

C 4 ^ {kl;^^c\d\ jD^){k^ - k\[d{pe\\ - Pex)/dpBp}^ The 
k‘^ term comes from the Hall term and is the stan¬ 
dard whistler wave, while the k‘j^ term is a whistler¬ 
like wave coming from the electron pressure anisotropy 
in generalized Ohm’s law. For eCGL, this is = 
(3/2)/3o/cywhere f3o = pojpB- Similarly, fol¬ 
lowing Ref. m, the intermediate frequency range gives 
a ;2 ^( 72 /^ 4 . In the k± ^ /c||, short wavelength, cold 
ion limit, this yields uu^ = {k‘^^k‘^c\d‘f/D‘^){(^dpe\\/dp)^ — 

[d(Pe±,Pe\\)/d(p,PB)]o}/[CA/D + {c\/D){dpeJ_/dpB)o + 
{dpe±/dp)o]. In the low Pq limit, the first term 
gives the standard kinetic Alfven wave, while the sec¬ 
ond is a kinetic Alfven-type wave arising from the 
pressure anisotropy. In eGGL, this wave has = 

{5/2)l3okjPd‘f{po/po)- 

There are many ways to test the dispersive wave model. 
For cold ions, dispersive waves from anisotropies per¬ 
sist. However, they vanish identically for cold electrons. 
The dispersive wave model predicts fast reconnection for 
eGGL but slow reconnection for iGGL, consistent with 
our simulations. Eq. @ implies there are dispersive 
waves without the Hall term when there is an equilib¬ 
rium pressure anisotropy, independent of the equations 
of state, consistent with previous studies [401156]. 

Interestingly, when Egedal’s equations of state [43l [44] 
are employed in simulations without the Hall term, re¬ 
connection is Sweet-Parker-like (J. Egedal, private com¬ 
munication). Thus, simply having an electron pressure 


anisotropy is insufficient to cause fast reconnection; the 
pressure anisotropy must have a particular form. Fluid 
modeling of other equations of state could provide insight 
about what physically sets the length of the current layer. 

We now show that electron pressure anisotropies can 
dominate the Hall term in real systems. First, we have 
performed particle-in-cell simulations with parameters 
similar to the fluid simulations, confirming that the CGL 
model (p|| oc p^/B^ and p± oc pB) is reasonably re¬ 
produced (plots not shown). Also, electron pressure 
anisotropies dominate the Hall term for the parameters 
of the simulations in Fig. Fig. [^f) shows a simulation 
of eCGL with the Hall term; the contribution to the re¬ 
connection electric field of the pressure anisotropy (solid 
line) dominates the Hall term (dashed line). 

We suspect electron pressure anisotropy dominates 
when dispersive wave terms due to the anisotropy dom¬ 
inate the standard whistler and kinetic Alfven waves in 
Eq. Q. In (72, the first term with d? gives the stan¬ 
dard kinetic Alfven wave. The second term with df is 
the most important term arising from the electron pres¬ 
sure anisotropy (by a factor of /d, which is small for many 
systems of interest). In the eg = e = 1 limit, a simple 
calculation reveals that the electron pressure contribu¬ 
tion of the kinetic Alfven wave is completely cancelled 
by part of the electron pressure anisotropy. This implies 
that it always dominates the Hall term when Tg > Ti 
with low Therefore, when > Ti, the anisotropy is 
the dominant mechanism for the entire parameter regime 
previously thought to be the kinetic Alfven regime of re¬ 
connection m - low [3, high in-plane based on Bj^qq, 
and strong guide field (but not strong enough to make the 
ion Larmor radius smaller than dg). Physically, Ti needs 
to be smaller than Tg because if it is large enough, it can 
dominate the electron pressure effect discussed here. 

There are physical systems where reconnection in this 
parameter regime could occur. The solar wind and some 
tokamaks are low {3 where significant guide fields are ex¬ 
pected and Tg > Ti is possible. 
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